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RATIONALE 


Even if the result of this meeting is an Ideal seeding device that gen- 
erates particles that exactly follow the flow and are of sufficient major source 
of error, I refer to particle counting bias wherein the probability of 
measuring velocity If It occurs Is a function of velocity* The error in the 
measured mean can be as much as 25% (ref. 1). 

Many schemes have been put forward to correct for this error, but there 
is not universal agreement as to the acceptability of any cone method. In par- 
ticular ft Is sometimes difficult to know if the assumptions required in the ana- 
lysis are fulfilled by any particular flow measurement system. 

In an effort to check various correction mechanisms In an Ideal way and to 
gain some Insight into how to correct with the fewest Initial assumptions, a 
computer simulation was constructed to simulate laser anemometer measurements in 
a turbulent flow. That simulator and the results of Its use are the topic of 
this paper. 


INTRODUCTION 

All measurements of mean quantities in a sparsely seeded turbulent flow 
using a laser anemometer generate a measured velocity probability function 
P m (v) that differs from the true Eulerian probability of interest P(v). The 
relation between the two is given by 


r (v) 

p n,< v > ' p <»> - 

m 

where r m (v) Is the measurement rate if 
the velocity is v. < > denotes 
average value. 


Any correction scheme's function is to eliminate the effect of r m (v). 
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SIMULATOR 

First a pseudo-continuous signal, the "hot wire” signal, is generated by 
passing a digital white noise signal through a digital filter* Shown here is a 
typical segment of the "hot wire” signal* 


Sample Hot Wire Velocity Signal 


Velocity 


The "laser anemometer” signal Is generated by e random sampling of the hot 
wire signal. The simulator is constructed so that the average measurement rate 
as a function of velocity can be set to be any desired function of the velocity. 
In this utudy the measurement probability, r m (v), ties set to be either a linear 
or quadratic function of the velocity magnitude. Shown here Is a block diagram 
of the particle sampling section. 


Flow Diagram of Particle Arrival Simulation 
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If one takes the ratio of the laser anemometer histogram to the hot wire 
histogram, the result should be the normalized rate corresponding to the velo- 
city (see the first figure). That ratio is shown here. Note the linear depen- 
dence of the rate on velocity. 







CORRECTION SCHEMES 

Three correction schemes from the literature were checked using the 
simulator. 

1) Sample and Hold. A continuous signal is generated by holding the last 
measured velocity until a new measurement arrives. That aew one Is the held... 
It has been predicted that In the limit of many measurements per flow correla- 
tion time, the continuous signal generates unbiased statistics. 

2) Periodic Measurement. Time is divided into intervals of constant 
length. If only the first measure in each interval is recorded, a periodically 
sampled signal is generated. It has been predicted that all bias vanishes when 
It is highly probable that there is a measurement in each interval. 

3) Dead Time. Any data recording device has a reset time during which no 
new measurement is recorded. If the dead time is small compared to the flow 
correlation time and if the particle rate is high enough, the bias should 
vanish. 

This figure shows each scheme graphically. 


Samplo Hot Wiro Velocity Signal 



Velocity 





The Dias in Che aean definitely does decrease as the mean rate per flow 
correlation time Increases. This Is shown In the next figure. 


Moan Measured Velocity ve. Sampling Rate 
V - 1.0. TI - 0.30 
Sample and Hold Analysis 

< v » 



Mean Sampling Rata Par Correlation Tima 


20 



If the aean particle arrival rate is held constant and if the periodic 
saapllng rate is varied, the bias in the aeasured mean does decrease but not to 
zero. See below. 


Measured Mean Velocity vs. Q 
V*. 997+/-. 008, TI". 300-*-/-. 003, santpl ing rate*l. 0 
Periodic Sampling Analysis 

< v » 



0 " ratio of aaMpllng period to corrolatlon tins 
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If the dead tiae Is kept smaller tian the flow correlation time, the bias 
Is supposed to decrease to cero as the particle arrival rate Increases. This Is 
borne out by the following figure. Although the figure doesn't show the bias 
going to tero, the measurements closely fit the theory that does go to zero. We 
haven't had enough computer time to check the high particle rate limit. 


Measured Mean Velocity ve. RT 
V * 1.0, TI * 0.30 
Dead Time Analysis 
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Recently Edwards and Meyers put forth a correction scheme that involved 
direct measurement of r m (v) (ref. 2). Simply, one measured that rate for each 
interval in the velocity histogram by counting the average number of measure- 
ments occurring in a small interval At after the appearance of a velocity in a 
given interval. It can be shown that the procedure's averaging result is exact 
in the limit of At going to zero. Unfortunately, it can also be shown that the 
relative measurement error goes to infinity as At goes to zero. The figure below 
shows the measured rate as a function of At. Note that as At increases, the 
measured rate tends toward the mean value, independent of which interval one 
starts in. 


Rato by Tim© of Arrival vs. Dolta - t 
Edwards Linear Method 


Rat. of Arrival 
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Edwards and Baratuci have invented a scheme that computes the limit as At 
goes to zero starting with relatively large values of At (ref. 3). A line is 
fit to the results for various At and the intercept is taken as the "correct" 
value. 


Rato of Particle Arrival vs. Oelta - t 
Velocity Histogram Bin § 5 
Edwards Method ( Linear) 

Rat* of Arrival 
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The following figure shows the measured rate versus velocity for a typical 
simulation with a quadratic dependence on velocity magnitude. The measured 
rates are very close to the values set by the simulation. 


Rate by Histogram vs. Velocity 
Edwards - Baratuci Method 


Rato by Histogram 



Weighted laoat 
S guars# Fit 


Vbar -1.0. T I - 0.30. Rt- - 1.0 


Velocity 


25 




i 

A 
. i 

,\ 

> » 



.# 


The following figure from reference 3 gives a summary of the results of 
our correction schemes for a set of simulations. McLaughlin and Tiederman (ref. 
1) and the histogram corrections are included only as a self consistency check 
of the simulator. The two Edwards correction schemes are here because In general 
one does not know the functional form of r m (v). Edwards (quadratic) is a general 
quadratic fit to r ffl (v) . 


Edwards - B.r.tuci Correction Method. 
Bate dependence on velocity is linear. 


Mean Turbulence 

Velocity Intensity 


Euler Ian 

0.998 

1 .008 

0.302 

£ .004 

Measured 

1.092 

£ .009 

0.269 

£ >004 

McLaughlin St 

1.000 

+ .012 

0.302 

£ .017 

Tiederman 

Hi stogram 

1.002 

t .010 

0.302 

£ .006 

( 1 i nssr ) 

Histogram 

1.003 

£ .012 

0.300 

£ .009 

(quadratic) 

Edwards 

1.000 

£ .027 

0.304 

£ .014 

( 1 i near ) 

Edwards 

1.010 

£ .023 

0.290 

+ .018 


(quadratic) 


Note : Vbar ** 1.0, TI - 0.30, Rf - 1.0. 

Note : Edwards (linear) means a linear lit was 

done on r (At,v ) and Edwards (quadratic) 
means a fit was done to a second order 
polynomial. The same explanation applies 
to the Histogram Methods. 

Note s Error bars are the standard deviation 

of the mean of the mean velocity obtained 
from twenty data ets of 2700 points each. 
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r * 


a, 






Euleri an 

1.000 ± .010 

0.301 £ .004 


Measured 

1.169 ± .013 

0.2S1 £ -304 

• 

McLaughlin It 

1.093 £ .013 

0.429 £ .01* 


Tiederman 




Histogram 

1.047 £ .029 - 

0.296 1 .030 


(linear) 




Histogram 

1.010 £ .019 

0.297 t *011 


(quadratic ) 




Edwards 

0.960 £ .097 

0.209 £ .033 

t 

(linear) 



j 

Edwards 

1.016 £ • 043 

0.274 £ .028 

<4 

(quadratic) 



A 

* 

* 

Note i Vbar * 

1.0, TI * 0.30, 

Rr c - 1.0, N * 2700. 

•; 

t 

» 

Note : Edwards 

(linear) means 

a linear fit was 

> » 

done on 

r (At,v.) and 
m k 

Edwards (quadratic) 



means a -fit was dons to a sscond order 
polynomial. Ths same explanation applies 
* ~ Histogram Methods. 

Not* : Error bars are ths ftandirj deviation 

a £ the mean a* the mee.i velocity obtained 
'from twenty data sets of 2700 points each. 
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| CONCLUSIONS 

i 

I 

| Sample and Hold eliminates bias for high 

j particle densities. 

1 

| Dead Time reduces bias for sampling rates tested 

1 C Prediction not checked due to large run time ] 

. Periodic Sampling with long periods reduces but 

does not eliminate velocity bias 

j 

New correction is excellent when the form of 
rate as a function of velocity is known 

Further work using "unknown" functional forms is 
i in process 
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